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We present the centrality dependent measurement of multiplicity and pseudorapidity distributions 
of charged particles and photons in Au + Au collisions at ^snn = 62.4 GeV. The charged particles 
and photons are measured in the pseudorapidity region 2.9 < rj < 3.9 and 2.3 < r/ < 3.7, respectively. 
We have studied the scaling of particle production with the number of participating nucleons and 
the number of binary collisions. The photon and charged particle production in the measured pseu- 
dorapidity range has been shown to be consistent with energy independent limiting fragmentation 
behavior. The photons are observed to follow a centrality independent limiting fragmentation be- 
havior while for the charged particles it is centrality dependent. We have carried out a comparative 
study of the pseudorapidity distributions of positively charged hadrons, negatively charged hadrons, 
photons, pions, net protons in nucleus-nucleus collisions and pseudorapidity distributions from p 
+ p collisions. From these comparisons we conclude that baryons in the inclusive charged particle 
distribution are responsible for the observed centrality dependence of limiting fragmentation. The 
mesons are found to follow an energy independent behavior of limiting fragmentation while the 
behavior of baryons seems to be energy dependent. 



PACS numbers: 25.75.Dw 

I. INTRODUCTION 

The STAR experiment Q at the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory has the unique capability of measuring 
charged particle and photon multiplicities, precisely 
and simultaneously, at forward rapidity. By using 
this capability we can carry out a systematic study 
of various aspects of charged particle and photon 
production in relativistic heavy ion collisions. 

The conventional way of describing particle pro- 
duction in heavy ion collisions is by measuring the 
particle density in pseudorapidity (rj). Within the 
framework of certain model assumptions, it pro- 
vides information on energy density, initial temper- 
ature and velocity of sound in the medium formed 
in the collisions 0- The widths of the pseudo- 
rapidity distributions are sensitive to longitudinal 
flow and re-scattering effects 0, The variation 
of particle density in r\ with collision centrality, ex- 
pressed in terms of the number of participating nu- 
cleons (iVpart) and/or the number of binary collisions 
(-Wcoii)j can shed light on the relative importance 
of soft versus hard processes in particle production. 
The particle density in pseudorapidity also provides 
a test ground for various particle production mod- 
els, such as those based on ideas of parton satura- 
tion p| and semi-classical QCD, also known as the 
color glass condensate (CGC) p. 

At RHIC, the particle production mechanism 
could be different in different regions of pseudorapid- 
ity. At midrapidity a significant increase in charged 
particle production normalized to the number of par- 
ticipating nucleons has been observed from periph- 
eral to central Au + Au collisions • This has been 
attributed to the onset of hard scattering processes, 
which scale with the number of binary collisions. 
However, the total charged particle multiplicity per 
participant pair, integrated over the whole pseudo- 
rapidity range, is independent of centrality in Au 
+ Au collisions |8|. In the framework of the color 



glass condensate picture of particle production [||, 
the centrality dependence of particle production at 
midrapidity reflects the increase of gluon density due 
to the decrease in the effective strong coupling con- 
stant. It will be interesting to see how the pho- 
ton and charged particle production scales with the 
number of participating nucleons and with the num- 
ber of binary collisions in a common r\ coverage at 
forward rapidity. The increase in particle produc- 
tion at midrapidity with increasing center of mass 
energy has been studied in detail at RHIC :8]. It is 
also of interest to see how particle production varies 
with center of mass energy at forward rapidity. The 
experimental data on hadron multiplicity and its en- 
ergy, centrality and rapidity dependence so far have 
been consistent with the approach based on ideas 
of parton saturation. Recently it has been argued 
that this onset of saturation occurs somewhere in 
the center of mass energy ( v / %n) range of 17 GeV 
to 130 GeV 9J. This is one of the reasons cited for 
having different mechanisms of particle production 
at RHIC and SPS. The present experimental data 
at y/sNN = 62.4 GeV may help to understand the 
transition energy for the onset of saturation effects 
in particle production. 

It has been observed that inclusive photon produc- 
tion (primarily from decay of 7r°) at y / SNN = 62.4 
GeV [l(| follows a centrality independent limiting 
fragmentation llj behavior. The inclusive charged 
particles at 19.6 GeV and 200 GeV have been ob- 
served to follow a centrality dependent behavior of 
limiting fragmentation §. It has been speculated 
that the baryons, an important constituent of in- 
clusive charged particles, are responsible for the ob- 
served difference between photons and charged par- 
ticles [lEil. The baryons coming from nuclear rem- 
nants and baryon transport, both of which change 
with centrality, may be the source of the central- 
ity dependent limiting fragmentation for inclusive 
charged particles. The role of a new mechanism of 
baryon production as discussed in Refs. [13, EH & l so 
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needs to be understood. A comparative study of 
limiting fragmentation of positively and negatively 
charged particles and photons at the same collision 
energy and pseudorapidity interval as provided by 
the present data will help to understand the sources 
responsible for the observed features. On the theo- 
retical side, reproducing the energy, centrality, and 
species dependence of limiting fragmentation ob- 
served in the experimental data can be a good test 
for various particle production models. One such at- 
tempt to explain the energy dependence of limiting 
fragmentation phenomena within the framework of 
CGC has been carried out in Ref. ^4|. The impor- 
tance of the limiting fragmentation curve on energy 
dependence of part icle production has been demon- 
strated in Ref. fig . 

Event-by-event measurements of photon and 
charged particle multiplicities can be used to study 
multiplicity fluctuations |l6j . Fluctuations in phys- 
ical observables in heavy ion collisions may provide 
important information regarding the formation of a 
Quark-Gluon Plasma and help to address the ques- 
tion of thermalization ^7|- The study of event-by- 
event fluctuations in the ratio of photon to charged 
particle multiplicities has also been proposed as a 
tool to search for production of Disoriented Chiral 
Condensates (DCCs) 0. 

In this paper we address some of the above physics 
issues through the first simultaneous measurement 
of the charged particle and photon multiplicities 
for Au + Au collisions at ^/snn = 62.4 GeV in 
the forward rapidity. The charged particles are de- 
tected using the Forward Time Projection Cham- 
ber (FTPC) and the photons are detected using the 
Photon Multiplicity Detector (PMD) in the STAR 
experiment |TL fl9L l20|| . 

The paper is organized as follows: In section II 
we briefly describe the detectors used for measuring 
the charged particle and photon multiplicities and 
the trigger detectors used for selecting the minimum 
bias data, used in the present analysis. In section 
III we present the details of data analysis from the 
FTPC and the PMD. In section IV we present the 
results in terms of multiplicity and pseudorapidity 
distributions of photons and charged particles, scal- 
ing of particle production with number of partici- 
pating nucleons and number of binary collisions and 
limiting fragmentation behavior for charged, neutral 
and identified particles. Finally we summarize our 
study in section V. 



II. DETECTORS 

The STAR experiment p] consists of several de- 
tectors to measure hadronic and electromagnetic ob- 
servables spanning a large region of the available 



phase space at RHIC. The detectors used in the 
present analysis are the FTPC, PMD, a set of trigger 
detectors used for obtaining the minimum bias data 
and the Time Projection Chamber (TPC), data of 
which is used to determine the collision centrality. 
The FTPCs, PMD, minimum bias trigger and colli- 
sion centrality selection are briefly described below. 



A. Forward Time Projection Chambers 

There are two FTPCs; they are located on each 
side of the nominal collision vertex, around the beam 
axis. They are cylindrical in structure with a diam- 
eter of 75 cm and 120 cm in length. Each FTPC has 
10 rows of readout pads, called padrows, which in 
turn are subdivided into 6 sectors with 160 pads per 
sector. The first padrow is located about 1.63 meters 
away on both sides from the center of the TPC (the 
nominal collision point). The sensitive medium is a 
gas mixture of Ar and CO2 in the ratio of 50%: 50% 
by weight. The FTPCs detect charged particles in 
the pseudorapidity region 2.5 < | 77 | < 4.0. The 
novel design of the FTPCs uses a radial drift field, 
perpendicular to the magnetic field, to achieve a two- 
track resolution up to 2 mm. This allows for track 
reconstruction in the environment of high particle 
density at forward rapidity. In the present analysis, 
the data from only the FTPC in the positive pseudo- 
rapidity region (2.9 < i] < 3.9) is used. Particle pro- 
duction models such as HIJING HU and AMPT H2 
show that about 6-7% of the total charged particles 
produced fall within the acceptance of each of the 
FTPCs. Further details of the design characteristics 
of the FTPC can be found in Ref. ]H|. 



B. Photon Multiplicity Detector 

The PMD is located 5.4 meters away from the cen- 
ter of the TPC (the nominal collision point) along 
the beam axis. It consists of two planes (charged 
particle veto and preshower) of an array of cellu- 
lar gas proportional counters. A lead plate of 3 ra- 
diation length thickness is placed between the two 
planes and is used as a photon converter. The sensi- 
tive medium is a gas mixture of Ar and CO2 in the 
ratio of 70%:30% by weight. There are 41,472 cells in 
each plane, which are placed inside 12 high voltage 
insulated and gas-tight chambers called supermod- 
ules. A photon traversing the converter produces an 
electromagnetic shower in the preshower plane, lead- 
ing to a larger signal spread over several cells as com- 
pared to a charged particle, which is essentially con- 
fined to one cell. The PMD detects photons in the 
pseudorapidity region 2.3 < r/ < 3.7. In the present 
analysis, only the data from the preshower plane has 
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TABLE I: Centrality selection, number of participating 
nucleons and number of binary collisions. 

~%CTOSS Section (iVpart) (iVcoll)" 



0- 


-5 


> 373 


347.3 


904.3 


5- 


10 


373-313 


293.3 


713.7 


10 


20 


313-222 


229.0 


511.8 


20 


30 


222-154 


162.0 


320.9 


30 


40 


154-102 


112.0 


193.5 


40 


50 


102-65 


74.2 


109.3 


50 


60 


65-38 


45.8 


56.6 


60 


70 


38-20 


25.9 


26.8 


70- 


-80 


20-9 


13.0 


11.2 



been used. From HIJING [2l| and AMPT we 
find that about 10-11% of the total photons pro- 
duced fall within the acceptance of the PMD. The 
photon production is dominated by photons from 
the decay of ir°s 0. HIJING calculations indicate 
that about 93-96% of photons are from inclusive 7T° 
decays. Further details of the design and character- 
istics of the PMD can be found in Ref . [2(j ■ 



C. Minimum bias trigger and collision 
centrality 

The minimum bias trigger is obtained using the 
charged particle hits from an array of scintillator 
slats arranged in a barrel, called the Central Trig- 
ger Barrel, surrounding the TPC, two zero degree 
hadronic calorimeters at ±18 m from the detector 
center along the beam line, and two Beam-Beam 
Counters |23J. The centrality determination in this 
analysis uses the uncorrected multiplicity of charged 
particles in the pseudorapidity region | 77 | < 0.5, as 
measured by the TPC |24|. Table [Q gives the per- 
centage cross section, the corresponding uncorrected 
multiplicity of charged particle tracks (iV2 ) in the 
pseudorapidity region | 77 | < 0.5, the number of par- 
ticipating nucleons (A^ par t) and the number of binary 
collisions (N co u) used in this paper. The number 
of participating nucleons and the number of binary 
collisions have been obtained from Glauber calcula- 
tions j24|. 



III. DATA RECONSTRUCTION 
A. Charged particle reconstruction 

The analysis of the data from the FTPC involves 
the following steps: (a) event selection, (b) pad- 
to-pad gain calibration, and (c) reconstruction of 
charged tracks. 

A total of 1.2 million minimum bias events, corre- 
sponding to 0-80% of the Au + Au hadronic interac- 



tion cross section, have been selected with a collision 
vertex position less than 30 cm from the center of the 
TPC along the beam axis. 

The calibration of the FTPC is done using a laser 
calibration system This system helps to cali- 

brate the drift velocity in the nonuniform radial drift 
field and also provides information for corrections to 
spatial distortions caused by mechanical or drift held 
imperfections. The localization of dead pads is done 
with pulsers and by an analysis of data to identify 
electronically noisy pads |25| . 

The reconstruction of experimental data involves 
two steps: (a) cluster-finding to calculate the track 
points from the charge distribution detected by the 
pads, (b) track-finding to group the track points 
of different padrows of the FTPC to form a track. 
The cluster-finding includes reading of the electronic 
signal data from the data acquisition system, look- 
ing for areas of nonzero charge (cluster), deconvo- 
lution of clusters, and then finding the point co- 
ordinates. This is followed by combining clusters 
from all padrows to form tracks using a suitable 
tracking algorithm pCSf. A track is considered valid 
if it consists of at least 5 found clusters and if its 
distance of closest approach to the primary ver- 
tex is less than 3 cm. The condition of having at 
least 5 found clusters for each track in the FTPC 
ensures a small contribution of split tracks. The 
split tracks contribution and background contamina- 
tion are primarily from 7 conversion electrons and 
positrons which are significantly reduced when we 
include those tracks which have transverse momen- 
tum in the range 0.1 < pt < 3 GeV/c in the analysis. 
The maximum percentage of split tracks was esti- 
mated from simulations to be ~1.5%. The relative 
amount of split tracks decreases as we go from from 
central to peripheral collisions. Two procedures are 
used to obtain the charged particle yields at all pt- 
The charged particle transverse momentum spectra 
are fitted by a power-law function in the range 0.1 
GeV/c < pt < 1 GeV/c and extrapolated to pt = 
GeV/c. The low pt yield is obtained from this 
extrapolation. The other procedure calculates the 
yield of charged particles for pt < 0.1 GeV/c by 
using the ratio of the yield in this range to to- 
tal yields from HIJING 0| simulations. Both these 
procedures resulted in similar correction factors of 
the order of 15% in the region 2.9 < rj < 3.9. 

The efficiency of charged particle reconstruction 
(e c h) as a function of pseudorapidity is estimated 
by embedding Monte Carlo charged tracks into 
real data and then following the full reconstruc- 
tion chain |27| . The reconstruction efficiency is 
obtained by dividing the number of reconstructed 
Monte Carlo tracks within an r\ bin by the total num- 
ber of embedded Monte Carlo tracks in the same ij 
bin. The charged particle reconstruction efficiencies 
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FIG. 1: (Color Online) Charged particle reconstruction 
efficiency (e c h) and purity of charged hadron sample (/ c h) 
in the FTPC as a function of pseudorapidity (77) for 
charged tracks with 0.1 GeV/c < pt < 3 GeV/c, for 
two centrality intervals. 



for central and peripheral collisions, for the 77 region 
studied, are shown in Fig. ^ The background con- 
tamination is obtained from detailed Monte Carlo 
simulation using the HIJING (version 1.382) event 
generator 12 ill and the detector simulation package 
GEANT |22[7 which incorporates the full STAR de- 
tector framework. The purity of the charged hadron 
sample (/ c h) in the FTPC for central and periph- 
eral collisions is also shown in Fig. ^ The errors on 
efficiency and purity values will be discussed later. 



B. Photon reconstruction 

The analysis of the data from the PMD involves 
the following steps: (a) event selection, (b) cell-to- 
cell gain calibration, and (c) reconstruction or ex- 
traction of photon multiplicity. 

A total of 0.3 million minimum bias events, corre- 
sponding to 0-80% of the Au + Au hadronic interac- 
tion cross section, have been selected with a collision 
vertex position less than 30 cm from the center of 
the TPC along the beam axis. The difference in the 
number of events for the PMD and FTPC analysis 
originates from the fact that, for the same period 
of data-taking, the PMD recorded fewer events and 
there was a need for a more stringent data clean-up 
procedure to remove events with pile-up-like effects. 

The cell-to-cell gain calibration was done by ob- 
taining the ADC distributions of isolated cells. The 
ADC distribution of an isolated cell may be treated 



as the response of the cell to charged particles |20J . 
For most of the cells this response follows a Landau 
distribution. We use the mean of the ADC distri- 
bution of isolated cells to estimate and correct the 
relative gains of all cells within each supcrmodule. 
The cell-to-cell gain variation is between 10-25% for 
different supermodules. 

The extraction of the photon multiplicity proceeds 
in two steps involving clustering of hits and photon- 
hadron discrimination. Hit clusters consist of con- 
tiguous cell signals. Photons are separated from 
charged particles using the following conditions: (a) 
the number of cells in a cluster is > 1, and (b) the 
cluster signal is larger than 3 times the average re- 
sponse of all isolated cells in a supermodule. The 
choice of the conditions is based on a detailed study 
of simulations |l0l |2Cj . The number of selected clus- 
ters, called 7 — like clusters (-/V 7 _iike), in different 
supermodules for the same 77 coverage are used to 
evaluate the effect of possible non-uniformity in the 
response of the detector. 

To estimate the number of photons (-/V 7 ) from the 



detected N~, 



-like 



clusters we evaluate the photon re- 



construction efficiency (e 7 ) and purity (/ p ) of the 



7 — like sample defined ;29| as e 7 = N^'^ /N 1 and 
/ P = N^l /N 7 -uke respectively. N]^ h is the num- 
ber of photon clusters above the photon-hadron dis- 
crimination condition. Both e 7 and / p are obtained 
from a detailed Monte Carlo simulation using HI- 
JING [2l| with default parameter settings and the 
detector simulation package GEANT |28| , which in- 
corporates the full STAR detector framework. Both 
e 7 and / p vary with pseudorapidity and centrality. 
This is due to variations in particle density, up- 
stream conversions and detector related effects. A 
photon should ideally create one cluster in the de- 
tector. However, it may give rise to more than one 
cluster (called split clusters) in the real experimental 
environment. These may happen because of conver- 
sions of the photon due to upstream materials in 
front of the PMD, or limitations of the clustering al- 
gorithm due to varying particle density. The highest 
occupancy of the PMD is about 12% and the max- 
imum percentage of split clusters is estimated to be 
9%. The photon reconstruction efficiency and the 
purity of the photon sample determined by means 
of simulations for central and peripheral collisions 
for the 77 region studied are shown in Fig. [21 The 
lower limit of photon px acceptance in the PMD is 
estimated from detector simulations to be 20 MeV/c. 



C. Systematic errors 

The systematic errors on the charged particle mul- 
tiplicity (-/V c h) are due to uncertainties in estimates 
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FIG. 2: (Color Online) Photon reconstruction efficiency 
(e 7 ) and purity of photon sample (/ p ) for PMD as a func- 
tion of pseudorapidity (rj), for two centrality intervals. 



of e c h and f^. The uncertainty in the estimates are 
obtained through simulations by varying the track 
quality cuts. The value of the maximal distance of 
closest approach of a track to the primary vertex is 
varied by 0.5 cm leading to a maximum error on N c h 
of ~6%. The minimum number of clusters to form 
a track was varied from 5 to 4. This led to an error 
on iV c h of The uncertainty in the correction 

factor to obtain the N c h yield for < 0.1 GeV/c is 
~8%. This also contributes to the total systematic 
errors. The total systematic error in iV c h is ~10% 
for all the centrality classes studied. The systematic 
error for the region rj > 3.6 is estimated to be about 
15%, due to larger uncertainity in the reconstruc- 
tion efficiency. This arises primarily due to uncer- 
tainity in realistic reproduction of electronic loss, at 
the extreme ends of the detector acceptance. This 
is estimated by studying the azimuthal dependence 
of charged particle density in a given rj window. 

The systematic errors on the photon multiplicity 
(Nj) are due to (a) uncertainty in estimates of e 7 
and / p values arising from splitting of clusters and 
the choice of photon-hadron discrimination thresh- 
old and (b) uncertainty in N 7 arising from the non- 
uniformity of the detector response primarily due to 
cell-to-cell gain variation. The error in iV 7 due to 
(a) is estimated from Monte Carlo simulations to 
be 9.8% and 7.7% in central and peripheral colli- 
sions, respectively. The error in JV 7 due to (b) is 
estimated using average gains for normalization and 
by studying the azimuthal dependence of the photon 
density of the detector in an 77 window to be 13.5% 
for central and 15% for peripheral collisions. The to- 
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FIG. 3: (Color Online) Minimum bias iV ch (2.9 < r/ < 
3.9 ) and iV 7 (2.3 < 77 < 3.7) distributions for Au + 
Au collisions at ^/sNN = 62.4 GeV. The charged particle 
and photon multiplicity distribution for top 5% central 
events are shown in open circles. The solid curve is the 
Gaussian fit to the data points. 



tal systematic error in iV 7 is ~17% for both central 
and peripheral collisions. 

The total errors on and iV 7 are obtained by 
adding respective systematic and statistical errors 
in quadrature and are shown in all the figures unless 
mentioned otherwise. The statistical errors are small 
and within the symbol sizes. 



IV. RESULTS AND DISCUSSION 
A. Multiplicity distributions 

The charged particle multiplicity (N c h) and pho- 
ton multiplicity (N^) are obtained event-by-event in 
the FTPC and the PMD following the analysis pro- 
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200 



FIG. 4: (Color Online) Correlation between average 
number of charged particles (N c h) and average number 
of photons (iV 7 ) within the common r\ range of FTPC 
and PMD 2.9 < r\ < 3.7 for different collision centrality 
classes in Au + Au collisions at ^/snn = 62.4 GeV. The 
solid line is a straight line fit to the data points (see text 
for details). 
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cedure described above. Fig. [31 shows the minimum 
bias distributions of iV c h and iV 7 for Au + Au colli- 
sions at ^/snn = 62.4 GeV. The distributions have 
a characteristic shape with a steep rise that corre- 
sponds to the most peripheral events. The plateaus 
in the photon and charged particle multiplicity dis- 
tributions correspond to mid-central events and the 
fall-off to the most central collision events. The 
shape of the curves in the fall-off region reflects 
the intrinsic fluctuations of the measured quantities 
and the limited acceptance of the detectors. The 
event-by-event charged particle and photon multi- 
plicity distributions for 0-5% central collisions are 
also shown. Gaussian fits to these distributions have 
been made. The values of the fit parameters for 
charged particles measured in 2.9 < r\ < 3.9 are: 
mean = 167 and a = 20; x 2 /ndf = 70.67/69. The 
values of the fit parameters for photons measured in 
2.3 <r/< 3.7 are: mean = 252 and a = 30; x 2 /ndf = 
37.3/34. The correlation between the average num- 
ber of charged particles and average number of pho- 
tons within the common pseudorapidity coverage of 
the FTPC and PMD (2.9 < 77 < 3.7) for different 
collision centrality classes in Au + Au collisions at 
•y/sNN = 62.4 GeV are shown in Fig. 21 The corre- 
lation between N&y and N 7 can be expressed as iV 7 
= (0.74±0.01)V ch - (3.57±0.83). This is shown as 
a straight line in the figure. The correlation reflects 
the variation of N 7 and N c h with collision centrality. 
The correlation coefficient is 1.01±0.01. 



FIG. 5: (Color Online) Variation of N c h normalized to 
the number of participating nucleon pair in the FTPC 
coverage (2.9 < r\ < 3.9) and AT 7 normalized to the num- 
ber of participating nucleon pair in the PMD acceptance 
(2.3 < rj < 3.7) as a function of -/V par t. The lower band 
shows the uncertainty in the ratio due to uncertainties 
in A^part calculations. 



B. Scaling of particle production 

After having discussed the event-by-event mea- 
surement of photon and charged particle multiplici- 
ties in the previous section, we now discuss the vari- 
ation of average (averaged over number of events) 
photon and charged particle multiplicities within the 
full coverage of the PMD and FTPC, respectively, 
with centrality. Collision centrality is expressed in 
terms of either number of participating nucleons or 
number of binary collisions. This will provide in- 
formation on the contribution of hard (pQCD jets) 
and soft processes to particle production at forward 
rapidity. The scaling of particle production with the 
number of participating nucleons indicates the dom- 
inance of soft processes while scaling with the num- 
ber of binary collisions indicates the onset of hard 
processes. At midrapidity the particle production at 
^/snn = 130 GeV and 200 GeV has been shown to 
scale with a combination of Af par t and A^oii [2J ■ Here 
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FIG. 6: (Color Online) Variation of N c ^ normalized to 
the number of collisions in the FTPC coverage (2.9 < 
r\ < 3.9) and iV 7 normalized to number of collisions, in 
the PMD coverage (2.3 < 77 < 3.7) as a function of iV co ii. 
The lower band shows the uncertainty in the ratio due 
to uncertainties in iV co ii calculations. 



we present the results on scaling of particle produc- 
tion at forward rapidity for Au + Au collisions at 
= 62.4 GeV. 

Figure [3] shows the variation of the total number 
of charged particles in the FTPC coverage (2.9 < 
77 < 3.9) and the total number of photons in the 
PMD acceptance (2.3 < 77 < 3.7), both normalized 
to -/Vpart, as a function of the collision centrality, ex- 
pressed by the number of participants. Higher iVp art 
values correspond to more central collisions, or col- 
lisions with smaller impact parameter. The charged 
particle yield per participating nucleon pair at for- 
ward rapidity decreases from peripheral to central 
collisions. The photon production per participant 
pair is found to be approximately constant with cen- 
trality in the forward 77 range studied. 

Figure shows the variation of the total number 
of charged particles normalized to the number of col- 
lisions in the FTPC coverage (2.9 < r; < 3.9) and the 
total number of photons normalized to the number 
of collisions in the PMD coverage (2.3 < 77 < 3.7) as 
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FIG. 7: (Color Online) dN/drj for charged particles and 
photons for Au + Au collisions at ^/ snn = 62.4 GeV for 
various event centrality classes. 



a function of the number of binary collisions. Higher 
^Vcoll values correspond to more central collisions, or 
collisions with smaller impact parameter. Both the 
charged particle yield and photon yield normalized 
to the number of binary collisions do not scale with 
the number of binary collisions at forward rapidity. 
The data value decreases from peripheral to central 
collisions. This indicates that the contribution of 
hard processes to particle production at forward ra- 
pidity is small. 



C. Pseudorapidity distributions 

So far we have discussed the multiplicities of pho- 
tons and charged particles over the full coverage of 
the detectors. In this section we study the varia- 
tion in particle density with 77. The results can then 
be directly compared to different models in order to 
understand the mechanism of particle production in 
heavy ion collisions at forward rapidity. 

Figure shows the pseudorapidity distributions 
of charged particles within 2.9 < 77 < 3.9 and pho- 
tons within 2.3 < 7/ < 3.7 for various event cen- 
trality classes. As expected the particle density in- 
creases with decrease in r\. Fig. |H1 shows the com- 
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FIG. 8: (Color Online) dN/dr/ for charged particles and 
photons for central and peripheral Au + Au collisions 
at -^/snn = 62.4 GeV compared to corresponding results 
from theoretical models. 



parison of pseudorapidity distributions for photons 
and charged particles for 0-5% and 40-50% central 
Au + Au collisions at v / snn = 62.4 GeV with the 
corresponding results from various theoretical mod- 
els. The HIJING model [2lJ is based on perturbative 
QCD processes which lead to multiple jet produc- 
tion and jet interactions in matter. HIJING seems 
to underpredict the measured photon multiplicity. 
However within the systematic errors it is difficult 
to make definitive conclusions. For charged parti- 
cles, HIJING fails to explain the r) distributions for 
central and peripheral collisions. The AMPT [22^ 
model is a multi-phase transport model which in- 
cludes both initial partonic and final hadronic inter- 
actions. For photons, the results from the AMPT 
model are in reasonable agreement with the data 
for central and peripheral events within the system- 
atic errors. For charged particles in central colli- 
sions, the results from AMPT explain the data at 
lower r\ and overpredict the charged particle yields at 
higher rj. The LEXUS |13 model is based on linear 
extrapolation of nucleon-nucleon collisions to high- 
energy nucleus-nucleus collisions. For charged par- 
ticles, the LEXUS model underpredicts the multi- 
plicity at lower r\ and agrees with experimental data 




FIG. 9: (Color Online) Pseudorapidity distributions of 
charged particles for various center of mass energies in 
Au + Au central collisions. The pseudorapidity distri- 
butions for v /i^ = 200 GeV, 130 GeV and 19.6 GeV 
are from the PHOBOS experiment ||. The solid lines 
are the results of the fits described in the text. 



at higher r\ for central collisions. It also underpre- 
dicts the charged particle yields for peripheral colli- 
sions. In summary, we observe that the photon mul- 
tiplicity within the systematic errors is reasonably 
well explained by HIJING and AMPT models. The 
detailed pseudorapidity dependence of the charged 
particle multiplicity is not reproduced by the above 
models. 

We have so far studied the pseudorapidity distri- 
bution of particles at forward rapidity at -^/snn = 
62.4 GeV. Now we will investigate the following: (a) 
the energy dependence of the shape of the r\ distribu- 
tion of charged particles available at various energies 
of Au + Au collisions in RHIC, and (b) try to esti- 
mate the full r) distribution for charged particles for 
■J snn = 62.4 GeV from the above study and com- 
pare to the present measurements. 

The full pseudorapidity distribution of charged 
particles at RHIC for central collisions can be 
parametrized by the following 3-parameter formula: 



dN_ 
dr\ 



C 



1 



exp^iw 



This formula is chosen to describe the central 
plateau and the fall off in the fragmentation region of 
the distribution by means of the parameters rjo and 
5 respectively. Using this formula we can describe 
the 200 GeV^ 130 GeV and 19.6 GeV pseudorapidity 
distributions of charged particles from the PHOBOS 
experiment [8j. The values of the parameters C, rjo 
and S are given in Table ITT1 and the fits to data are 
shown in Fig. EI The value of rjo is found to increase 
with increasing a/snn- The value of the parameter <5 
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TABLE II: Parameters C, 770 and 8 for different 



C 



'7o 



19.6 382 ± 33 2.16 ± 0.17 0.7 ± 0.06 

62.4 (interpolated) 458 ± 40 3.08 ± 0.35 0.69 ± 0.06 
130 580 ± 21 3.59 ± 0.076 0.66 ± 0.05 

200 667 ± 22 3.80 ± 0.082 0.71 ± 0.06 
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FIG. 10: (Color Online) Half width at half maximum of 
the pseudorapidity distributions (r/h) of charged particles 
as a function of total charged particle multiplicity (Nt) 
normalized to the center of mass energy. The Au + Au 
collision data are from the PHOBOS 81 experiment and 
p + p collision data are from the ISR |31f experiments. 



is found to be independent of energy within errors. 
This feature is another way of testing the concept 
of limiting fragmentation, which will be discussed 
later. Using the average value of S and interpolating 
the value of 770 to 62.4 GeV we are able to predict 
the full pseudorapidity distribution for charged par- 
ticles at 62.4 GeV. This is shown as solid curve in 
Fig. together with our measured charged particle 
data for 62.4 GeV at forward rapidity. The dashed 
curves represent the error in obtaining the full pseu- 
dorapidity distribution for charged particles using 
the interpolation method described. 

We have also studied the widths of the pseudo- 
rapidity distributions of charged particles at RHIC 
and compared them to those from p + p collisions 
at ISR [3l| ■ In Fig. ^] we show the variation of 
the half width at half maximum (7^) of the charged 



particle pseudorapidity distributions as a function 
of total charged particle multiplicity normalized to 
the center of mass energy (JVt / -^/snn) for P + P and 
Au + Au collisions. The data shown is for various 
centrality classes in Au + Au collisions || and for 
various intervals of observed total multiplicity in p + 
p collisions. We observe that the half width at half- 
maximum obeys an interesting scaling law in p + p 
collisions and is found to depend on a single vari- 
able (JVt / ^/snn) ■ In Au + Au collisions this scaling 
seems to be valid for 200 GeV and 130 GeV. Al- 
though the width decreases with A't/v'snn for 19.6 
GeV, the data lies below the higher energy data un- 
like the energy independent behavior observed in p 
+ p collisions. This may reflect the change in the 
mechanism of particle production over the full pseu- 
dorapidity range as we increase the ^/snn from 19.6 
GeV to -^/snn > 130 GeV in Au + Au collisions at 
RHIC. 



D. Energy dependence of particle multiplicity 

The energy dependence of charged particle yields 
at midrapidity has been studied at RHIC §]. Here 
we present the results on the energy dependence of 
particle yields at forward rapidity and compare them 
with yields at midrapidity. 

Figurc lTTl shows the charged particle pseudorapid- 
ity distribution scaled by the number of participat- 
ing nucleon pairs at midrapidity (77 = 0) and for- 
ward rapidity (77 = 3.0) as a function of ^/snn for 
central collisions at RHIC. The data for charged par- 
ticles at ysiii = 19.6 GeV, 56 GeV, 130 GeV and 
200 GeV at 77 = 3.0 are from the PHOBOS and 
BRAHMS [331 experiments. The data for charged 
particles at midrapidity are the averages of the val- 
ues from the 4 RHIC experiments. The charged par- 
ticle production at 77 = 0, can be expressed as 



dN/dr] 



0.5N, 



part 



= 1.75(±0.25) + 0.017(±0.005) In [v^nn] 
-0.00003(±0.00002)(ln[ % /iNN]) 2 . 



The charged particle production at rj = 3.0, can be 
expressed as 



dN/dr] 
0.5JV part 



-0.03(±0.13) + 0.028(±0.004) In [v^nn] 
-0.00007(±0.00002)(ln[ v /i^]) 2 . 



The ratio of charged particle production at 77 = 
to that at 77 = 3.0 decreases from a factor 4 to 1.3 
as ^/snn increases from 19.6 GeV to 200 GeV. The 
photon result at ^/snn = 62.4 GeV for 77 = 3.0 is 
also shown. The photon yields at other ^/snn val- 
ues at forward rapidity and midrapidity are not yet 
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FIG. 11: (Color Online) dN/dr) per participating nu- 
cleoli pair at midrapidity (77 = 0) and forward rapidity 
(rj = 3.0) for various center of mass energies for central 
collisions. The data for charged particles at ^/snn = 19.6 
GeV, 56 GeV, 130 GeV and 200 GeV at r, = 3.0 are from 
the PHOBOS [| and BRAHMS [13 experiments. The 
data for charged particles at midrapidity are the aver- 
ages of the values from the 4 RHIC experiments. The 
photon yield at ^/ snn = 62.4 GeV is also plotted. For 
comparison the results from a model based on parton 
saturation expected at high-density QCD |33| are also 
shown. The solid lines are polynomial fits to the values 
from the QCD model. There is no prediction for a/snn 
= 62.4 GeV available from this model. 



available at RHIC. The photon production at ^/snn 
= 62.4 GeV is about 35% lower than the charged 
particle production for the same energy at 77 = 3.0. 
The charged particle yield at 77 = 3.0 for ^/snn = 
62.4 GeV is a factor 1.6 and 1.9 lower compared to 
the corresponding yields at 130 GeV and 200 GeV 
and a factor 3.0 higher than the charged particle 
yields at 19.6 GeV. For comparison, also shown in 
Fig. ^2 are the results from a model based on par- 
ton saturation, which is expected in high-density 
QCD 33]. The results from the model agree with 
the measured charged particle yields at midrapidity 
for all energies at RHIC. However, the model's pre- 
diction for forward rapidity at the lowest energy (22 
GeV) is lower compared to data (19.6 GeV). There 
is no prediction for ^/snn = 62.4 GeV available from 
this model. It would be interesting to have the pre- 
dictions to understand the transition energy for the 
onset of saturation effects at RHIC. 
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FIG. 12: (Color Online) Ratio of iV c h to iV 7 for 0-5% and 
40-50% central Au + Au collisions at ^/s^ = 62.4 GeV 
as a function of 77. Results from HIJING are also shown 
for comparison. The lower band reflects the common 
errors in ratio for the two centrality classes. 



E. Comparison of N c h and N-, 

The STAR experiment at RHIC has the unique ca- 
pability to study the yields of charged particle and 
photons at forward rapidity. Fig. El shows the ra- 
tio of N ch to 7V 7 for 0-5% and 40-50% central Au 
+ Au collisions at ^/snn = 62.4 GeV as a function 
of 77 in the common 77 coverage of the FTPC and 
the PMD. The ratio is around 1.4 for central colli- 
sions and 1.6 for peripheral collisions within 3.0 < 
77 < 3.6. The results from HIJING indicate simi- 
lar values. The correlated systematic errors, mainly 
arising due to uncertainties in the Monte Carlo de- 
termination of reconstruction efficiencies and nor- 
malization errors, are not plotted on the data points 
and are shown as a shaded band. The photon pro- 
duction is dominated by photons from the decay of 
7r°s The charged particle yields have a sub- 

stantial contribution from baryons at forward rapid- 
ity |34j . Apart from the kinematics, this may be the 
reason for higher charged particle yields compared 
to photons. In the future, event-by-event study of 
7V c h and N 7 correlations in common 77 and <p cover- 
age of the FTPC and the PMD can be used to look 
for possible formation of disoriented chiral conden- 
sates |Tg| . 



F. Energy dependence of limiting 
fragmentation 

Continuing our discussion on particle density in 
77, we now present results on the longitudinal seal- 



13 



g.2 



LO 

o 



12 



lo 
o 







=A 0-5% 130 GeV (BRAHMS) 
P 0-6% 130 GeV (PHOBOS) 
|0 0-6% 19.6 GeV (PHOBOS) 



0-5% 62.4 GeV (STAR) 
0-6% 200 GeV (PHOBOS) 



%0 




O 0-5% 62.4 GeV (STAR) 
□ 0-5% 17.3 GeV (WA98) 
O 0-5% 19.6 GeV (WA93) 



-2.5 



-1.5 -1 

^beam 



-0.5 



FIG. 13: (Color Online) The top panel shows the varia- 
tion of dNch/dr; normalized to N pax t with r\ - ybeam for 
different collision energies for central collisions. The bot- 
tom panel shows the variation of dN-y/dr] normalized to 
A?part with r\ - ybeam for different collision energies for 
central collisions. 



ing of particle production in heavy ion collisions. It 
has been observed that the number of charged par- 
ticles produced per participant pair as a function of 
V - ybeam, where ybeam is the beam rapidity, is inde- 
pendent of beam energy 0, llfjfl . This phenomenon 
is known as limiting fragmentation 11]. Here we 
present the results on energy dependence of limiting 
fragmentation at 62.4 GeV for charged particles and 
photons produced in Au + Au collisions. In the sub- 
sequent sections we discuss the centrality and species 
dependence of this scaling. 

In Fig. E| we present the energy dependence of 
limiting fragmentation for inclusive charged particles 
and photons. The charged particle pseudorapidity 
distribution for central (0-5%) Au + Au collisions 
at yj snn = 62.4 GeV is compared to the charged 
particle pseudorapidity distributions from PHOBOS 
for central (0-6%) collisions at 19.6 GeV, 130 GeV 
and 200 GeV @ and charged particle pseudorapidity 
distribution from BRAHMS for central (0-5%) col- 
lisions at 130 GeV The photon pseudorapidity 
distributions for central (0-5%) Au + Au collisions 
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FIG. 14: (Color Online) The top panel shows the vari- 
ation of dTVch/dr? normalized to 7V pa rt with r\ - ybeam 
for the central and peripheral collisions. Bottom panel 
shows the variation of dN^/dr/ normalized to iVpart with 
T) - ybeam for the central and peripheral collisions. Also 
shown are the charged particle and photon yields in p + 
p and p + p collisions. 



at ^/snn = 62.4 GeV is compared with central (0- 
5%) photon data for Pb + Pb collisions at 17.3 GeV 
from the WA98 experiment [2!j and 19.6 GeV cen- 
tral (0-5%) S + Au collision data from the WA93 
experiment [3^|. We observe in Fig.^|that the SPS 
and RHIC (62.4 GeV) photon results are consistent 
with each other, suggesting that photon production 
follows an energy independent limiting fragmenta- 
tion behavior. The charged particles at 62.4 GeV 
also show an energy independent limiting fragmen- 
tation behavior. 



G. Centrality dependence of limiting 
fragmentation 

Recently there have been contradictory results re- 
ported from inclusive charged particle measurements 
regarding the centrality dependence of the limit- 
ing fragmentation behavior. Results from PHO- 
BOS show a centrality dependence @, while those 
from BRAHMS show a centrality independent be- 
havior [33. Here we present the results on the 
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centrality dependence of limiting fragmentation for 
charged particles and photons at snn = 62.4 GeV. 

In Fig. E| we show the centrality dependence 
of limiting fragmentation for charged particles and 
photons. The charged particle pseudorapidity dis- 
tributions for 0-5% and 40-50% central Au + Au 
collisions at -^/snn = 62.4 GeV have been compared. 
We observe, at forward rapidity, the charged par- 
ticle yield normalized to the number of participat- 
ing nucleons as a function of rj - ybeam is higher 
for peripheral collisions compared to central colli- 
sions, whereas within the measured rj range of 2.3 
to 3.7, the photon yield normalized to the number 
of participating nucleons as a function of rj - ybeam 
is found to be independent of centrality. The de- 
pendence of limiting fragmentation on the collision 
system is most clearly seen in the comparison be- 
tween results from heavy ion collisions with those 
from p + p and p + p collisions 3j| . We observe in 
Fig. I14l that the photon results in the forward rapid- 
ity region from pp collisions at ^/snn = 540 GeV are 
in close agreement with the measured photon yield 
in Au + Au collisions at v / snn= 62.4 GeV. However 
the p + p and p + p inclusive charged particle results 
are very different from those for Au + Au collisions 
at v / snn= 62.4 GeV. It may be mentioned that the 
photon yield is dominated by photons from decay 
of 7r°s 10]. The presented photon results and their 
comparison with nucleon-nucleon collisions indicate 
that in the rj region studied, there is apparently a 
significant charged baryon contribution in nucleus- 
nucleus collisions. Similar centrality dependent be- 
havior of limiting fragmentation for charged parti- 
cles was also observed by PHOBOS |S|]. The central- 
ity dependence of limiting fragmentation in charged 
particles has been speculated to be due to nuclear 
remnants and baryon stopping 0, Il2| . The central- 
ity independent limiting fragmentation for photons 
has been attributed to mesons being the dominant 
source of photon production ^(j ■ HIJING calcula- 
tions indicate that about 93-96% of the photons are 
from 7T° decays. 

In order to understand the role of nuclear rem- 
nants and baryon stopping in the observed central- 
ity dependent behavior of limiting fragmentation 
of charged particles, we have studied the limiting 
fragmentation for positively and negatively charged 
hadrons separately. The contribution from protons 
coming from beam remnants can be understood by 
studying the limiting fragmentation of positively 
charged hadrons. In Fig. ED we have plotted 
normalized to the number of participating nucleons 
for 40-50% and for 0-5% central collisions for pos- 
itively (h + ) and negatively charged (h~) hadrons. 
In addition to the systematic errors discussed ear- 
lier, and shown in the figure, there is an error due 
to the uncertainty in the charge determination. The 
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FIG. 15: (Color Online) Variation of dN c h/drj normal- 
ized to JVpart with rj - ybeam for central and peripheral 
collisions for positively charged hadrons (h + ) and nega- 
tively charged hadrons (h~). 



uncertainty has been studied by embedding charged 
Monte Carlo tracks into real data and then follow- 
ing the full reconstruction chain. This error was ob- 
tained as a function of rj. It is defined as the ra- 
tio of the total number of embedded charged tracks 
whose charge has been reconstructed incorrectly, to 
the total number of charged tracks embedded. The 
error in charge determination was found to increase 
from 2% at rj = 2.9 to 15% at T] = 3.9. We find 
that both h + and h~ show a centrality dependent 
limiting fragmentation behavior. When compared 
to the centrality independent limiting fragmenta- 
tion behavior for photons (Fig. I14f> and to results 
from nucleon-nucleon collisions (Fig. 1 1 4|) , our mea- 
surements indicate that baryon transport at forward 
rapidity also plays an important role in the observed 
centrality dependent behavior of limiting fragmenta- 
tion for charged particles. We find that the ratio for 
yields of h + from peripheral to central collisions in- 
creases from 1.17±0.06 at rj = 3.0 to 1.61±0.07 at 
rj = 3.8 (closer to beam rapidity). The values for 
h~ are 1.16±0.06 at rj = 3.0 and 1.51±0.07 at rj = 
3.8. From these values we find that the increase in 
the ratio with rj seems to be somewhat weaker for 
h~ compared to h + . However, within the system- 
atic errors, it is difficult to conclude on the role of 
the beam remnants (beam protons in h + ) in the cen- 
trality dependent behavior of limiting fragmentation 
for charged particles at forward rapidity. 
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FIG. 16: (Color Online) Variation of dN^/dr] normal- 
ized to Afp ar t with r\ - ybeam from AMPT model [2^| cal- 
culations for various center of mass energies in central 
collisions (top panel) and central and peripheral colli- 
sions at x /snn= 19.6 GeV (bottom panel). 



Energy and centrality dependence of limiting frag- 
mentation for charged particles can be a test for par- 
ticle production models. We have observed that par- 
ticle production models such as HIJING and AMPT 
are not able to describe fully the r] distribution of 
charged particles at forward rapidity. However, it is 
interesting to investigate whether they can qualita- 
tively reproduce the limiting fragmentation features 
of experimental data. Our calculations show that 
in the HIJING and AMPT models the charged par- 
ticles show energy independent limiting fragmenta- 
tion. The centrality dependent behavior of limiting 
fragmentation for charged particles is more clearly 
observed in the AMPT model than in HIJING. In 
Fig. we only show the results from the AMPT 
model. These results are for the -^/snn values of 
19.6 GeV, 62.4 GeV and 130 GeV Au + Au collisions 
at 0-3 fm and 9-12 fm impact parameter. For the 
centrality dependence we only show the results for 
V s nn= 19.6 GeV, the energy at which the centrality 
dependent effect is most prominent in the data 
The AMPT model has qualitative limiting fragmen- 
tation features similar to those of experimental data 
(shown in Fig. I14|) . We find in the model that the 



central yields, when normalized to number of partic- 
pating nucleons, are also lower than the correspond- 
ing peripheral yields at forward rapidity when -q is 
shifted by the beam rapidity. 



H. Identified particle limiting fragmentation 

The observation of centrality dependent and en- 
ergy independent limiting fragmentation for inclu- 
sive charged particles, along with the centrality and 
energy independent limiting fragmentation for pho- 
tons (presented in previous sections), motivates us 
to study the limiting fragmentation of identified par- 
ticles. 

The top panel in Fig. El shows the charged pion 
rapidity density in central Au + Au collisions at 
RHIC |37j, Pb + Pb collisions at the SPS [H and 
Au + Au collisions at AGS Q. Also shown is the 
estimated 7r° rapidity density from the present mea- 
surement of the photon rapidity densit y a t ^/snn = 
62.4 GeV, all as a function of y-ybeam [Hi- We ob- 
tained the ratio of the photon to tt° yields from HI- 
JING. This ratio is used to estimate the 7r° yield 
from the measured photon yield. The results indi- 
cate that the pion production in heavy ion collisions 
in the fragmentation region agrees with the energy 
independent limiting fragmentation picture. 

The bottom panel of Fig. 1 171 shows the net proton 
(p - p) rapidity density in central Au + Au colli- 
sions at RHIC 34] energies and Pb + Pb collisions 
at SPS energies. For AGS energies pi lioT| we 
plot only the proton rapidity density in Au + Au 
collisions. Since the anti-proton yields are very low 
(p/p ~ 2 x 10~ 4 at top AGS energy), the proton 
rapidity density reflects the net proton rapidity dis- 
tribution. The net protons violate the energy depen- 
dence of limiting fragmentation. These results show 
that baryons and mesons differ in the energy de- 
pendence of limiting fragmentation. The results for 
identified particles, along with the centrality depen- 
dence of limiting fragmentation for inclusive charged 
hadrons, and the centrality independence of limiting 
fragmentation for identified mesons, shows that the 
baryon transport in heavy ion collisions plays an im- 
portant role in particle production at forward rapid- 
ity. The results also show that although baryon stop- 
ping is different in different collision systems, the pi- 
ons produced at forward rapidity are not affected by 
baryon transport. The limiting fragmentation study 
for net protons may also indicate the validity of a 
baryon junction picture If the- baryon numbers 
are carried by the valence quarks, then at forward 
rapidity the baryons should also follow an energy in- 
dependent limiting fragmentation behavior, like pi- 
ons (originating from valence quarks). This may 
indicate that the baryon number is not carried by 
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In summary, we have presented charged particle 
and photon multiplicity measurements at RHIC 
in the pseudorapidity regions 2.9 < r) < 3.9 and 
2.3 < rj < 3.7, respectively. The pseudorapidity 
distributions of charged particles and photons for 
Au + Au collisions at ^/snn= 62.4 GeV have 
been obtained for various centrality classes and 
compared to results from different models. Charged 
particle and photon production normalized to the 
number of participating nucleon pairs and to the 
number of binary collisions has been studied. The 
photon multiplicity, within the systematic errors, 
seems to scale with the number of participating 
nucleons, while the charged particle multiplicity 
does not. Both the photon and charged particle 
production at forward rapidity do not scale with 
number of binary collisions. This indicates that 
the particle production at forward rapidity is not 
dominated by a contribution from hard processes. 
Charged particle and photon distributions at ^/snn 
= 62.4 GeV are both observed to be consistent 
with the energy independent limiting fragmentation 
scenario. The photon production is observed to 
follow a centrality independent limiting fragmenta- 
tion scenario, while the charged particles follow a 
centrality dependent behavior. Comparison of the 
pseudorapidity distributions of positively charged 
particles, negatively charged particles, pions, and 
distributions from p + p collisions, indicate that 
the baryons are responsible for the centrality de- 
pendent limiting fragmentation behavior of charged 
particles. The study of limiting fragmentation for 
pions and net protons show that mesons follow 
energy independent limiting fragmentation, whereas 
baryons do not. 



FIG. 17: (Color Online) The top panel shows the vari- 
ation of pion rapidity density normalized to N pa ,n with 
y _ ybeam for central collisions at various collision en- 
ergies. Also shown is the estimated dN^o /dy obtained 
from dN 1 /dy normalized to A part . The bottom panel 
shows the variation of net proton rapidity density nor- 
malized to iVpart with y - yboam for central collisions at 
various collision energies. 



the valence quark, which is suggested in the baryon 
junction picture, where the baryon number resides 
in a non-perturbative configuration of gluon fields, 
rather than in the valence quarks. 
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